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ABSTRACT 

We report the C0( J =1 — 0) observations of the Whirlpool Galaxy M51 using both the Combined 
Array for Research in Millimeter Astronomy (CARMA) and the Nobeyama 45m telescope (NR045). 
We describe a procedure for the combination of interferometer and single-dish data. In particular, we 
discuss (1) the joint imaging and deconvolution of heterogeneous data, (2) the weighting scheme based 
on the root-mean-square (RMS) noise in the maps, (3) the sensitivity and uv-coverage requirements, 
and (4) the flux recovery of a combined map. We generate visibilities from the single-dish map and 
calculate the noise of each visibility based on the RMS noise. Our weighting scheme, though it is 
applied to discrete visibilities in this paper, should be applicable to grids in ww-space, and this scheme 
may advance in future software development. For a realistic amount of observing time, the sensitivities 
of the NR045 and CARMA visibility data sets are best matched by using the single dish baselines only 
up to 4-6 kX (about 1/4-1/3 of the dish diameter). The synthesized beam size is determined to conserve 
the flux between synthesized beam and convolution beam. The superior ww-coverage provided by the 
combination of CARMA long baseline data with 15 antennas and NR045 short spacing data results in 
the high image fidelity, which is evidenced by the excellent overlap between even the faint CO emission 
and dust lanes in an optical Hubble Space Telescope image and PAH emission in an Spitzer Sfxm image. 
The total molecular gas masses of NGC 5194 and 5195 {d = 8.2 Mpc) are4.9xl09MQ and T.SxIO^Mq, 
respectively, assuming the CO-to-H2 conversion factor of Xqo = 1-8 x 10^*'cm~^[K • km/s]~^. The 
presented images are an indication of the millimeter- wave images that will become standard in the next 
decade with CARMA and NR045, and the Atacama Large Millimeter/Submillimeter Array (ALMA). 

Subject headings: techniques: interferometric — techniques: image processing — galaxies: individual 
(NGC 5194, NGC5195, M51) 
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1. INTRODUCTION 

Interferometers have an intrinsic limitation, namely, 
the problem of missing information. An interferome- 
ter records the target Fourier components of the spa- 
tial emission distribution, but an interferometer with a 
small number of antennas (N) can collect only a lim- 
ited number, N{N — l)/2, of the Fourier components 
instantaneously. In addition, the finite diameter of each 
antenna limits the minimum separation between anten- 
nas, which, in turn, imposes a maximum size on an ob- 
ject that the interferometer can detect. The zero-spacing 
data (i.e. zero antenna separation data) carry the impor- 
tant information of the total flux, and this information 
is always missing. The incomplete Fourier coverage (uv- 
coverage) also degrades the quality of image. Deconvo- 
lution schemes have been developed to extrapolate the 
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observed uv data to estimate the missing information, 
however, the performance is poor for objects with high 
contrast, such as spiral arms and the interarm regions of 
galaxies. 

The small- iV problem is particularly severe in millime- 
ter astronomy, though it is greatly reduced with the 
15-element Combined Array for Research in Millime- 
ter Astronomy (CARMA). CARMA combines the previ- 
ously independent Owens Valley Millimeter Observatory 
(OVRO) array {N — 6) and Berkeley-Illinois-Maryland 
Association (BIMA) array {N = 10 - reduced to 9 for 
CARMA). The number of antenna pairs, or baselines, is 
increased to 105 from the previous values of 15 (OVRO) 
and 45 (BIMA), providing a substantial improvement in 
uv coverage. In most observatories, a few array configu- 
rations are used to increase the number of baselines. The 
uv coverage from one CARMA configuration is equiva- 
lent to that from seven configurations with a 6-element 
array. CARMA ensures the unprecedented uv coverage 
in millimeter interferometry compared to previous mm- 
wave arrays. 

Single-dish telescopes complement the central uv cov- 
erage and provide short baselines, including the zero- 
spacing baseline. The combination of interferometer and 
single-dish data is not trivial, though several methods 
have been suggested. Existing methods can be cate- 
gorized into three types. The first method produces 
visibilities from a single-dish map (IVogel et al.l [1981 
iTakakuwal I2003D and adds single-dish and interferome- 
ter da ta in the uv domain. Pety & Rodrfguez-Fernandei 
(j2010t) discussed a mathematica formalism. One is- 
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sue faced when utilizing this method, the diiJicuhies 
of which are discussed in ([Heifer et aLl 120031) . has been 
the weighting of the two set s of data in com bina- 
tion. Rodrig uez-Fernan dez. Pe tv fc GuethI (|2008l ) and 
iKurono. Morita fc Kamazakil 12009') manually set the 
single-dish weight relative to the weight of the interfer- 
ometer to improve the shape of the synthesized beam. 
In this paper, we suggest a new weighting scheme based 
solely on the quality of the single-dish data. In our 
method, the single-dish weight is independent of the in- 
terferometer data and is intrinsic to the single-dish ob- 
servations. It naturally down-weights (up-weights) the 
single-dish data when its quality is poor (high). In 
the appendix, we discuss the sensitivity matching which 
makes the combination most effective. 

The second type of combination method co-adds two 
sets of data in the image domain fSt animirovic et al.l 
Il999t ) . This approach produces a joint dirty image and 
synthesized beam by adding the single-dish map and 
the interferometer dirty image, and single-dish beam and 
interferometer synthesized beam, respectively. The joint 
dirty image is then deconvolved with the joint dirty 
beam. This technique was ad opted for the BIM A-SONG 
survey (iHelfer et al.l [20031). iCornwelll ([T988I) . and re- 
cently iStanimirovic et al.l ([20021 ) . also discussed a non- 
linear combination technique through joint deconvolu- 
tion with the maximum entropy method ( MEM). 

The third method was introduced by iWeifi et al. I 
((20010 and operates in the Fourier plane. The de- 
convolved interferometer map and single-dish map are 
Fourier transformed and then the central uu-space from 
interferometer data is replaced with single-dish data. 

This paper describes the observations, data reduc- 
tion, and combination of CARMA and NR045 data of 
M51. Our procedure unifies the imaging techniques for 
interferometer mosaic data, heterogeneous array data, 
and the combined data of single-dish and interferome- 
ter. E arlier data reducti on and results have been pub- 
lished ([Koda et al.ll2009t ). The method and results are 
the same, but we have re-calibrated and reduced the 
entire data set using higher accuracy calibration data. 
In and g3l we describe the CARMA and NR045 ob- 
servations and calibration. The deconvolution (such as 
CLEAN) is detailed in fJD for three cases: (a) homoge- 
neous array, single-pointing observations, (b) heteroge- 
neous array, single-pointing observations, and (c) het- 
erogeneous array, mosaic observations. The weighting 
scheme in co-adding the images from a heterogeneous ar- 
ray (with multiple primary beams) is discussed in §4.41 
The result from this subsection is also essential for the 
combination of interferometer and single-dish data. The 
conversion of a single-dish map to visibilities is explored 
in 33 and discusses the resultant map and image fi- 
delity. A summary of the requirements of single-dish ob- 
servations for the combination are explained in ^j7| Com- 
ments on other combination methods are given in J}8l 
The summary is in ^ and sensitivity matching between 
single-dish and interferometer observations is discussed 
in Appendix [Cl 

2. CARMA 

® The synthesized beam is the instrumental point spread function 
for the aperture synthesis array; also know as the "dirty beam". 



2.1. Observations 

High resolution observations of the Whirlpool galaxy 
M51 in the CO(J = 1 — 0) line were performed with 
the Combined Array for Research in Millimeter Astron- 
omy (CARMA) during the commissioning phase and in 
the early science phase during the CARMA construc- 
tion (2006-2007). CARMA is a recently-developed in- 
terferometer, combining the six 10-meter antennas of 
the Owens Valley Radio Observatory (OVRO) millime- 
ter interferometer and the nine 6-meter antennas of the 
Bcrkcley-Illinois-Maryland Association (BIMA) interfer- 
ometer. The increase to 105 baselines provides superior 
uu-coverage and produces high image fidelity. The C 
and D array configurations are used. The baseline length 
spans over 30-350 m (C array) and 11-150 m (D array). 

The observations started with the heterodyne SIS re- 
ceivers from OVRO and BIMA. The typical system tem- 
perature of these original receivers was ~ 200 K in 
double-side band. The receivers of the 15 antennas were 
being upgraded one antenna at a time during the period 
of observations, but the process was not completed be- 
fore these observations finished. The system temperature 
of the new replacement receivers is typically ^ 100 K at 
115 GHz. 

The first-generation CARMA digital correlators were 
used as a spectrometer. They had three dual bands (i.e. 
lower and upper side bands) for all 105 baselines. Each 
band had five configurations of bandwidth - 500, 62, 31, 
8, and 2 MHz - which have 15, 63, 63, 63, and 63 chan- 
nels, respectively. We switched the configuration of band 
1, 2, 3 between (band 1, 2, 3) = (500, 500, 500) for gain 
calibration quasar observations and (band 1, 2, 3) — (62, 
62, 62) for target integrations. This "hybrid" configura- 
tion ensures both a sufficient detection of the gain cali- 
brator 1153-1-495 with the total 3 GHz bandwidth (i.e. 3 
bands x 2 side bands x 500MHz bandwidth) and a suffi- 
ciently wide velocity coverage for the main galaxy NGC 
5194. The total bandwidth is 149.41 MHz after drop- 
ping edge 6 channels at each side, which could be noisier 
than the central channels. The companion galaxy NGC 
5195 was not included in the velocity coverage, although 
it was detected in the NR045 map ( ^13. II) . 

The hybrid mode observations require a special calibra- 
tion for amplitude and phase offsets between bands and 
between configurations. We observed a bright quasar by 
changing the correlator configurations in time sequence: 
1. (band 1, 2, 3) = (500, 500, 500), 2. (62, 62, 62), 3. 
(500, 62, 62), 4. (62, 500, 62), and 5. (62, 62, 500). Each 
configuration spends 5 min on integration, and the whole 
sequence takes 25min integration in total. We used the 
bright quasars 3C273, 2C279, or 3C345, depending on 
availability during the observations. For any pair of band 
and bandwidth, this sequence has simultaneous integra- 
tions which can be used to calibrate the phase offset and 
amplitude scale between bands. The calibration observa- 
tions took typically 45 min including the radio pointing 
and antenna slew. These integrations were used for pass- 
band calibration as well. 

An individual observation consisted of a 4-10 h track. 
The total observing time (after flagging tracks under bad 
weather) is about 230 h (^ 30 tracks). A typical track 
starts with radio pointing observations of a bright quasar 
available at the time, then observes a fiux calibrator (e.g. 
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a planet), and repeats the 25 min observing cycle of gain 
calibrator (~ 5 min) and target (20 min including an- 
tenna slew for mosaic). The passband/hybrid observa- 
tions were performed at the middle of a track when M51 
is at a high elevation 80 deg). At such a high eleva- 
tion, each antenna slew between M51 and the calibrator 
takes a considerable amount of time. Observing a pass- 
band calibrator at a lower elevation avoids this loss. The 
system temperature (Tsys) was measured every gain cal- 
ibrator cycle, and the atmospheric gain variation is cor- 
rected real-time using Tsys. We observed 1153+495 as a 
gain calibrator. 

The telescope pointings were corrected every 4 h during 
the night and every 2 h during daytime. The last ^ 10 
tracks of the 30 total tracks also included an ad ditional 
optica l pointing procedure developed by Cord er et al.l 
(|2010f l. The optical procedure can operate during day- 
time, as well as at night, and a pointing correction was 
made every gain calibration cycle. This method measures 
the offset between radio and optical pointing vectors at 
the beginning of track (which is stable over periods much 
longer than the typical observation). During the ob- 
serving cycle of gain calibrator and target, the pointing 
drift, typically several arcsec per hour, is adjusted using 
a bright star close to the gain calibrator using an optical 
camera. The overhead of the pointing adjustment is less 
than 1 min. 

We mosaiced the entire 6'.0 x 8'. 4 disk of M51, with the 
disk defined by optical images and shown in Figure [U in 
151 pointings with Nyquist sampling of the 10m antenna 
beam (FWHM of 1 arcmin for the 115GHz CO J=l-0 
line). Ideally, every pointing position would be observed 
every M51 observing cycle (~ 20 min duration) to main- 
tain a uniform data quality and uv coverage across the 
mosaiced area. However, the overhead for slewing is sig- 
nificant for the large mosaic. It is as long as 6 sec per 
slew, and about 15 min total for 151 pointings. We there- 
fore observed every third pointing (total ~ 50 pointings) 
in each observation cycle to reduce the overhead. Three 
consecutive cycles cover all 151 pointings. Each track 
started from a pointing randomly chosen from the table 
of the 151 pointings, which helps the uniform data qual- 
ity among pointings. The resultant CARMA uv coverage 
is very similar at all pointings, and an example of the uv 
coverage at the central pointing is in Figure [2j 

The primary flux calibrators, Uranus, Neptune, and 
MWC349, were observed in most tracks. We monitored 
the flux of gain calibrator 1153+495 every month over the 
course of the observations. The flux of 1153+495 varied 
slowly between 0.7 and 1.3 Jy. The CARMA observatory 
is separately monitoring the flux variations of common 
passband calibrators, and our flux measurements are con- 
sistent with the observatory values. 

2.2. Calibration 

The data were reduced and calibrated using the 
Multichannel Image Reconstruction, Image Anal- 
ysis, and Display (M IRIAD) software package 
(|Sault. Teuben fc Wrightl 1^95). We developed ad- 
ditional commands/tasks to investigate and to reduce 
the large amount of data effectively, and to combine 
interferometer and single-dish data. 

The initial set of calibrations are the required routines 
for most CARMA data reduction. First, we fiag the data 



with problems such as antenna shadowing and bad Tsys 
measurements. Second, we apply the correction for vari- 
ation of optical fiber cable length, namely line length cor- 
rection. CARMA is a heterogeneous array of two types 
of antennas (i.e., 6m and 10m), and the optical fiber ca- 
bles that connect the antennas to the control building 
are mounted differently for the 10m and 6m dishes. The 
time variations of the cable lengths due to thermal ex- 
pansion are therefore different, which results in phase 
wraps in the baselines between 6m and 10m antennas. 
The changes of the cable lengths were monitored to an 
accuracy of 0.1 pico-second by sending signals from the 
control building and measuring their round-trip travel 
time. The changes are stored in MIRAD data and are 
used for the line-length correction. Third, we smooth the 
spectra with the Hanning window function to reduce the 
high side-lobes in raw spectra from the digital correla- 
tors. The spectral resolution is lowered by a factor of 2 
and becomes 1.954 MHz (5.08 km/s at the C0( J = 1-0) 
frequency) . 

Calibrations for passband and hybrid correlator config- 
uration were made using the sequence of hybrid configu- 
ration observations described in §2.11 We first separate 
500 MHz and 62 MHz integrations from the sequence and 
make two MIRIAD data sets containing only 500MHz 
data or 62 MHz data. These data sets are used to derive 
and apply passbands. The passband calibration removes 
the phase and amplitude offsets among Band 1, 2, and 
3 in the 500 and 62 MHz modes. An offset/passband 
calibrator is significantly detected even in 10 sec integra- 
tion both in the 62 MHz and in the 500 MHz mode. We 
derive the phase offset and amplitude scale between the 
500 MHz and 62 MHz modes by comparing the visibili- 
ties from the two modes on the 10 sec integration basis, 
and averaging them over time to derive single values for 
the phase offset and amplitude scale. We applied these 
calibrations to the entire track, which removes the phase 
and amplitude offsets between gain calibrator and target 
integrations. Errors of the hybrid calibration are small 
compared to the other errors and are only a few percent 
in amplitude and a few degrees in phase. 

The last set of calibrations includes the standard phase 
calibrations to compensate for atmospheric and instru- 
mental phase drifts. We did not use the gain calibra- 
tor integrations with large phase scatters (due to bad 
weather) and flagged the target integrations in the cy- 
cles immediately before and after the bad gain data. 
The absolute fluxes of the gain calibrator were measured 
monthly against a planet ( )j2.ip and were applied to tar- 
get data. 

The resulting la noise level of the CARMA data is 27 
mjy/beam in each lOkms""'^ channel. 

3. NOBEYAMA RADIO OBSERVATORY 45M TELESCOPE 

3.1. Observations 

We obtained total power and short spacing data 
with the 5 x5-Beam Array Receiver System (BEARS; 
iSunada et~a l. 2000) on the Nobeyama Radio Observa- 
tory 45m telescope (NR045). The FWHM of the NR045 
beam is 1 5'^ at 115 GHz. W e configured the digital spec- 
trometer ()Sorai et al.l I2OOO0 to 512 MHz bandwidth at 
500 kHz channel resolution. This is wide enough to cover 
the entire M51 system (both NGC 5194 and 5195). Han- 
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ning smoothing was applied to reduce the side-lobe in 
channel, and therefore, the resolution of raw data is 1 
MHz. 

We scanned M51 in the RA and DEC directions 
using the On-The-Fly (OT F) map pi ng technique 
(iMangum. Emerson fc GreisenI 120071 : iSawada et alJ 
|2008). We integrated OFF positions around the galaxy 
before and after each ~ 1 min OTF scan. A scan starts 
from an emission-free position at one side of the galaxy 
and ends on another emission-free position at the other 
side. Spectra are read-out every 0.1 second interval 
during the scan. The receiver array was rotated by 7 deg 
with respect to the scan directions, so that the 25 beams 
draw a regular stripe with a 5" separation. In combining 
the RA and DEC scans, the raw data form a lattice with 
5" spacing. This fine sampling, with respect to the beam 
size of 15", is necessary in reproducing the uv data up to 
the 45m baseline (i.e. the diameter of NR045), since we 
need the Nyquist sampling (5.96") of \go/D = 11.92", 
where Aco is the wavelength (= 2.6 mm) and D is th e 
antenna diameter (jMangum. Emerson &: GreisenI [20071 ) . 
If the sampling is coarser than 5.96", the aliasing effect 
in the Fourier space contaminates even shorter baseline 
data significantly. For example, if the sampling spacing 
is only 10.3" (i.e., typical sam pling in past NR045 
observations; iKuno et al]|2007[ ). the uv data down to 
the ^ 7 m baseline is contaminated (Figure [3]), and 
cannot be combined with the interferometer data. 

The typical system temperature in double side band 
was ~ 320 K. The pointing of the telescope was checked 
every ~ 45 min and was accurate to within ^ 2-3". 
BEARS is an array of double-side band (DSB) receivers 
and provides the antenna temperature (DSB) in DSB. 
The upper/lower side band ratio, namely the scaling fac- 
tor, was measured by observing Orion IRC2 using both 
BEARS and the single-side band (SSB) receiver SlOO and 
taking the ratios of the two measurements. The error in 
the measurements is a few percent. The total observing 
time under good weather conditions is about 50 hours. 

3.2. Calibration 

The ON/OFF calibration to account for the sky back- 
ground level was applied after the observations. We in- 
terpolated between two OFF-sky integrations before and 
after each OTF scan (^ 1 minute long), which reduced 
non-linear swells in the spectral baselines significantly. 
We used the NOSTAR data reduction package devel- 
oped at the Nobeyama Radio Observatory (jSawada et al.l 
[200l . converted the flux scale from T*(DSB) of BEARS 
to T*(SSB) of SlOO, subtracted linear spectral baselines, 
and flagged bad integrations. 

The 5" lattice of data from the observations was re- 
gridded with a spheroidal smoothing function, resulting 
in a final resolution of 19.7". We used the grid size of 
5.96", which is the Nyquist sampling of the 45 m spacing 
in Fourier space; this pixel scale is necessary to prevent 
artifacts from the aliasing effect (i j3.ip . 

We made maps of the RA and DEC scans separately. 
The two maps were co-added after subtracting spatial 
baselines in each scan direction to reduce systematic er- 
rors in the scan direction. Note that for OTF mapping, 
the sharing of an OFF among many ON scans may in- 
troduce noise correlations, primarily at small spatial fre- 
quencies in the Fourier space. lEmerson fc Grave. (,1988.) 



reduced such correlated noise using the basket-weave 
method, which down-weights the data at small spatial 
frequencies in the scan directions when the RA and DEC 
maps are added. We compared the spatial-baseline sub- 
traction and basket-weave methods, and found that both 
diminish the large-scale noise well. The difference was 
subtle, but the former gave a slightly smaller RMS noise, 
and thus, we decided to use the spatial-baseline method. 

The antenna temperature T^(SSB) was converted to 
the main beam temperature Tmb, using the main beam 
efficiency of r;„,b = 0.4 and r„,b = T;(SSB)/r;„,b. 

The flux of the flnal NR045 map is consistent with 
most previous measurements within a typical error of 
millimeter-wave measurements (10-20%). It is com- 
pared with four other results: an image from the Na- 
tional Radi o Astronomy Obs ervatory 12 m telescope 
(NRA012; 'Heifer et al. "2003^, two previous measure- 
ment s at NR045 (.Nakai et al. .19941 : [Matsushita et al.l 
mig) and our new CARMA data (^. The fluxes 
from iHelfer et al.] (|2003D . [Matsushita et al.l (|1999[ ). and 
the new CARMA observations are 94%, 95%, and 93% 
of that of the new NR045 map, respectively. For the 
comparisons, we re-sampled the new map to match the 
area coverage of the other maps. For the comparison 
with CARMA, the CARMA wii-distribution is ge nerated 
from the new NR045 map (as discussed in ij5.11 but for 
Hatcreek, OVRO, and CARMA primary beams), and the 
positive fluxes (above about Act) in the dirty maps are 
compared to measure the flux ratio. We used a Gaus- 
sian taper (FWHM=20") to make the dirty maps, which 
roughly reproduces the weight distribution of the NR045 
data. 

Only the map of iNakai et aTl ()1994l . distributed 
through iKuno et al.l ()200 7*)) shows a significant discrep- 
ancy: a factor of 1.82 higher total ffux than the new 
NR045 map. We attribute this discrepancy to an error 
in the old map, since all other measurements are consis- 
tent. Among these measurements, we decided to rely on 
the CARMA flux because we had the best understand- 
ing of the process of flux calibration, and because it is 
based on multiple flux calibrations over the duration of 
the observations. We scaled the ffux of the NR045 map 
to match the CARMA ffux (i.e., multiplied 0.93). 

The la noise level of NR045 data is 14.7 mK in 
rQ(SSB), 36.7 mK in Tmb, and 155 mJy/beam in 
lOkms^^ channel. 

4. IMAGING HETEROGENEOUS-ARRAY MOSAIC DATA 

We use MIRIAD for joint-deconvolution of multi- 
pointing CARMA and NR045 data. The method and 
algorithm for mosaic data with a homogencovis array are 
described in Sault. Stavcly-Smith & Brouw (1996). Our 
imaging involves two additional complications: a hetero- 
geneous array, and combinations with single-dish data, 
as well as mosaicing. We describe the essence of joint- 
deconvolution using MIRIAD, with an emphasis on the 
case of CARMA and NR045. 

Two points are of particular importance: the treat- 
ment of different primary beam patterns, and the weights 
of the data from the different primary beam patterns and 
from the single-dish. Here, we illustrate these two points, 
and define our notations. 

Correction for primary beam attenuation is simple for a 
homogeneous array. All antennas have the same primary 
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beam pattern P{l,m), and the primary beam correction 



I{l,m) = 



I{l,m) 
P{l,m)' 



(1) 



where the primary-beam corrected irnage is denoted / 
and the uncorrected image is denoted /. _The sky coordi- 
nates are (I, m). The uncorrected image / has two advan- 
tages: the synthesized beam B (i.e., point spread func- 
tion, PSF) and noise level are position- invariant, which 
simplifies the process of deconvolution 

For a heterogeneous array, the differences between pri- 
mary beam patterns have to be taken into account. For 
example, CARMA has three baseline types (i.e., antenna 
pairs), which result in three primary beam patterns - 
called "H" for Hatcreek (6m-6m dish pair), "O" for 
OVRO (lOm-lOm), and "C" for CARMA baseline types 
( 6m- 1 0m). Using appropriate weights Wh, Wq, and Wc 
(23), the images from "O", "H", and "C" baselines can 
be added as 



-TH -ro Pc 

= WkIh + WoIo + WcIc- 



(2) 



The weight is a function of position {l,m). The co- 
added image has been corrected for primary beam atten- 
uation. In the co-added plane, the synthesized beam pat- 
tern B and noise level are position-variant, which com- 
plicates the deconvolution. 

4.1. Homogeneous Array, Single- Pointing Data 

Traditionally, the imaging of interferometer data has 
been performed as follows. A set of one dirty map 
I'^™(l^m) and one synthesized beam pattern B{l,m) is 
made from visibilities. The dirty map I'^™ is deconvolved 
with B. For example, the deconvolution scheme CLEAN 
replaces the pattern B{1 — lo,m — mo), centered at an 
emission peak at {Iq, mo), with an ellipsoidal Gaussian 
to reduce the sidelobes of B. CLE AN_ usually runs in 
the /'*™ domain; the synthesized beam B and noise level 
a are position-invariant and their treatments are simple. 
The CLEANed image /™p is corrected for primary beam 
attenuation (cq. ([T])), providing the final map /™p. We 
note again that primary beam uncorrected and corrected 
images (of any kind) are differentiated with "bar" (e.g., 

jdm jdm jmp jmp-)_ 

The deconvolution is also possible in the J'^™ domain. 
The synthesized beam pattern B and noise level are not 
position-invariant. Thus, we define a position-variant 
synthesized beam pattern. 



B{l,m\ lQ,mo) 



B{1 — lo,m — mo) 
P(/,m) '' 



(3) 



centered at {lo,rnQ), and a position-variant noise level 
a/ P{l,m). Emission peaks are searched on a basis of 
signal-to-noise ratio. In MIRIAD, a set of primary-beam 
corrected / and uncorrected B is calculated from visibil- 
ities, and the command "mossdi" (i.e., CLEAN) calcu- 
lates B with eq. ^ at peak position (Zo,mo). 

4.2. Heterogeneous Array, Single Pointing Data 



The deconvolution in the image domain is applicable 
to heterogeneous array data. The joint dirty image I'^™ 
is defined as a linear summation of three dirty maps 
/^™, /g™, and /g™ (eq. The corresponding syn- 

thesized beam B is also a linear summation with the 
same weights. 



B(l,m; lQ,mo) = Wiiil,m] 



Woil,m: 



Bh{1 -lo,m- 


nio) 


Pnihrn) 




Boil -lo,m- 


mo) 


Po{l,m) 




Bc{l -lQ,m- 


mo) 



(4) 



+ Wc{l,m, 

Pc{l,mj 

The MIRIAD command "invert" with the mosaic op- 
tion outputs a set of one joint dirty map /'^™ (primary 
beam corrected) and three synthesized beams Bn, Bq, 
and Be (uncorrected). The command "mossdi" finds a 
peak emission in J'^™, and calculates B at its position 
with eq. (g]). 

In the case of a heterogeneous array, such as CARMA, 
the primary beam correction always needs to be applied 
to the dirty map. Thus, even for single-pointing obser- 
vations, we always use "options=mosaic" for "invert". 

4.3. Heterogeneous Array, Mosaic Data 

The deconvolution of mosaic data with a heterogeneous 
array is a further extension of the same procedure. Eq. 
(21) is extended as 



I{1, m)^Y. ^".f ^ = E ^^-'P^^'P 



(5) 



where the summation is taken for all baseline types h and 
pointings p. Wt.p is a weight for b and p. In practice, P is 
truncated at some radius, and only a subset of pointings 
contribute to a given position. 

The joint synthesized beam is defined as in eq. (|4]), 
but includes all pointings. In the case of the CARMA 
M51 observations, the command "invert" with "op- 
tion^mosaic" outputs one joint dirty map and 453 syn- 
thesized beams (= 3 baseline types x 151 pointings). A 
joint synthesized beam B is calculated with the 453 syn- 
thesized beams for every emission peak in I*^™. 

The spatial resolution is calculated by taking a 
weighted average of all 453 synthesized beams using Wb 
(b =11, O, C) and by filing a Gaussian. In theory, the 
sizes of the synthesized beams are different among the 
pointings, since the uv coverage is not exactly the same 
for all of the pointings. In practice, we designed the ob- 
servations to provide uniform uv coverage for all point- 
ings ( ii2.ip . We therefore adopt a single beam size over 
the whole mosaic. 

4.4. Weighting 

The noise level is position-dependent, a / P{1, m), in the 
image domain. Therefore, the weights W are defined as 

2 



Wh{l, m) (X 



Pb{l,m) 



(6) 



for 6 = H, O, and C, and are normalized as Wu -I- Wq + 
Wc = 1 at each position {l,m). The theoretical noise 
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a depends on baseline type h, and is the same as the 
imaging sensitivity AS*' discussed below. 

4.4.1. Thermal Noise and Its Coejficient 

Two sensitivities, fringe sensitivity AS'* and 
imaging sensitivity AS' [= cr], are important 
(jTavlor. Carilli fc PerlevI 119991 see their section 9). 
The fringe sensitivity is a sensitivity per visibility. The 
theoretical sensitivity for each visibility is calculated 
with the system temperature Tgys, bandwidth B, and 
integration time of the visibility tvis as 

where 



The aperture efficiency r]a and collecting area A of an- 
tennas i and j have a relation with the beam solid angle 
f^A given by l/{riaA) = fl^/X^. The Boltzman constant 
is fee and the last term l/^/2riq is due to the backend 
(i.e., dig itizer and correlator), a nd rjq is the quantum ef- 
ficiency (jRohlfs fc Wilsoiill200CII ). C^- is approximated as 
a constant for a homogeneous array, since the parameters 
are very similar for all antennas. In the case of a hetero- 
geneous array, Cy depends on baseline type. Parameters 
are listed in Table [TJ 

The imaging sensitivity is a root-mcan-square (RMS) 
noise in a final image, and depends on control parameters 
(see Appendix VK\ e.g., natural and uniform weighting). 
If the natural weighting is employed, the imaging sensi- 
tivity is simply a statistical summation of fringe sensitiv- 
ities, 1/(AS")^ — Efcl/(AS'*)^. For a homogeneous array 
{Cij is constant), it is 

AS'Na]=C,,Y^^&, (9) 

assuming that Tgys is a constant during observations. 
The total integration time is ftot = -^visivis, where A^vis 
is the number of visibilities. 

5. THE COMBINATION OF NR045 WITH CARMA 

The NR045 image is converted to visibilities and com- 
bined with CARMA data in uv space. Here we dis- 
cuss four steps for combination: (1) generating visibili- 
ties from the single-dish image, (2) the calculation of the 
weights of the single-dish visibilities, in the same form as 
interferometer visibilities, (3) the determination of syn- 
thesized beam size, and (4) an imaging/deconvolution 
scheme. A flow chart of the procedure is in Figure |4l 

5.1. Converting the NR045 Map to Visibilities 

To produce NR045 visibilities, we first deconvolve a 
NR045 map with a NR045 point spread function (PSF), 
multiply a dummy primary beam, generate a Gaussian 
visibility distribution, and calculate the amplitude and 
phase of the visibilities from the deconvolved, primary- 
beam applied NR045 map (Figure |4]). The following 
sections describe these steps. 

One limitation arises from the current software, though 
it should be easily modified in future software develop- 
ment. NR045 visibilities must have the same form as 



those of interferometers, and therefore, a dummy pri- 
mary beam needs to be applied to the NR045 map. 

5.1.1. Deconvolution with the NR045 Beam 

A NR045 map is a convolution of a true emission 
distribution with a point spread function (PSF). In the 
case of OTF mapping (i jS.ll) . the PSF is not literally the 
NR045 beam, but is a convolution of the NR045 beam 
and the spheroidal function which is used to re-grid the 
observed data to a map grid (i i3.2|) . The intrinsic Gaus- 
sian FWHM of the NR045 beam is 15", and is degraded 
to 19.7" after the re-gridding. The NR045 map needs to 
be de-convolved with this PSF. 

Figure [5] shows the sensitivity (noise) as a function of 
uu-distance (baseline length). It has a dependence on the 
Fourier-transformed PSF FT{PSF} as cx l/v/FT{PSF} 
(see Appendix [C]). The standard deviation of FT{PSF} 
is o-F = 3.9 kA for a Gaussian PSF with the FWHM of 
19.7". Thus, the noise increases significantly beyond 4-6 
kA (i.e. -\/2crF). Figure [5] shows that the NR045 sensitiv- 
ity is comparable to that of CARMA up to 4-6 kA and 
deviates beyond that. With the resultant sensitivities, 
we decided to fiag the data at > 4kA. The long baselines 
have negligible effects if we us e only the w eight based 
on sensitivity (i.e. robust=+2, lBriggsl[T995D . but could 
introduce an elevated error when robust < +2. 

CARMA and NR045 are complementary in terms 
of ttf-coverage and se nsitivity (Figure [2] and [5]). 
iKurono. Morita &: Kamaz aki (2009) suggested that the 
single-dish diameter should be 1.7 times as large as the 
minimum baseline of interferometer data, which is ~ 18 
meters in our case. However, we seem to need a 45m class 
telescope to satisfy the sensitivity requirement within re- 
alistic observing time. The sensitivity matching between 
NR045 and CARMA data is discussed in Appendix [C] 

5.1.2. Applying a Dummy Primary Beam 

The imaging tasks in MIRIAD assume that all visi- 
bilities are from interferometric observations, and apply 
a primary beam correction in the process of imaging. 
Consequently, the NR045 visibilities need to be attenu- 
ated by a pseudo primary beam pattern Pn- The choice 
of Pn is arbitrary, and we employ a Gaussian primary 
beam with the FWHM of 2 arcmin. Pn is multiplied to 
the deconvolved NR045 map at each of the 151 CARMA 
pointings separately. Since the map will be divided by 
Pn during the deconvolution, the choice of Pn does not 
affect the result. However, it is safer to use Pn at least 
twice as large as the separation of the pointings, so that 
the entire field is covered at the Nyquist sampling (or 
over-sampling) rate. 

We note that this multiplication of a primary beam in 
the image domain is equivalent to a convolution in the 
Fourier domain. It smoothes the sensitivity distribution 
in uw-space, and therefore, the weight discussed in ^5.21 
The size of the primary beam in Fourier space is only 
1/6 of that of the NR045 beam. Therefore, this effect 
should be small and negligible. 

5.1.3. Generating a Gaussian Visibility Distribution 

The distribution of visibilities in uv space should repro- 
duce the NR045 beam (more precisely, the PSF in i)5.ip 
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as a synthesized beam in image space. The Fourier trans- 
formation of a Gaussian PSF is a Gaussian. Therefore, 
visibihties are distributed to produce a Gaussian density 
profile in uv space. The size of the Gaussian distribution 
is set to reproduce the beam size of 19.7". We manually 
add a visibility at (u, v) = (0,0), so that the zero-spacing 
is always included. The number of visibilities TVvis and 
integration time per visibility iint are control parameters, 
and are discussed in ^5: 



5.1.4. Resampling 

Prom the Gaussian visibility distribution and the pri- 
mary beam attenuated maps, the visibility amplitudes 
and phases are derived, which gives the NR045 visibili- 
ties. 

5.2. Theoretical Noise and Other Parameters 

The relative weights of the CARMA and NR045 visi- 
bilities are important for proper combination. MIRIAD 
requires a weight (sensitivity) per individual visibility 
for imaging, and we calculate the weight based on the 
RMS noise of a NR045 map. For the interferometer 
data (i j4.4.1l) . we start from the fringe sensitivity and 
calculate the imaging sensitivity 5' by summing up the 
S^s of all visibilities. Here, we start from the RMS noise 
of a map (i.e., 5') and determine and its coefficient. 

The theoretical noise of a single-dish map, in main 
beam temperature Tmb, is 



AT, 



mb 



(10) 



where r]q and rj^h are the quantum efficiency of the spec- 
trometer and the main beam efficiency of the antenna, 
respectively. B and ttot are the bandwidth and total inte- 
gration time, respectively. [Note that the contribution to 
the noise from the OFF posi tion integrat ions should be 
negligible in OTF mapping (jSawada et a l. 2008)]. The 
total integration time (per point) of the NR045 map is 
derived from the RMS noise in the map using this equa- 
tion. 

The imaging sensitivity, corresponding to eq. (jO]), is 
calculated by converting the unit of eq. pH)) from Kelvin 
to Jy, 



AS*' 



2kv 



sys 



Comparing with eq. ©, we obtain 

2fcB 1 

~ J — • 

VmhVaA riq 

The fringe sensitivity per visibility should be 



sys 



(11) 



(12) 



(13) 



where the integration time per visibility is tvis — ^tot/A'vis 
and iVvis is the number of visibilities. The ^vis value 
should be set (arbitrary) to a small number, so that A'vis 
becomes large enough to fill the uv space. We set ^vis — 
0.01 sec and iVvis = 42075. 

Conceptually, we can understand the meaning of the 
NR045 visibilities by comparing the definition of fringe 
sensitivities (eqs. [7l and [T3| . They are the ones observed 



virtually with two identical NR045 antennas. The two 
antennas can physically overlap (in our virtual observa- 
tions), so they can provide uv coverage down to zero 
spacing. The beam shape of the NR045 dish plays the 
role of synthesized beam, but not primary beam. The 
primary beam shape is arbitrarily defined by Pn ~ if we 
seek a meaning, it corresponds to the beam shape of 
small patches within the NR045 dishes. 

There is one caveat when this weighting method is ap- 
plied with the current version of MIRIAD. MIRIAD is 
designed for an array with the same backend for all visi- 
bilities, and therefore, it neglects the l/^/2rjq term from 
Cij (eq. [S]). It defines an alternative parameter. 



JYPERK 



2k^ 



(14) 



which is stored in data header. The weights (AS*^) are 
calculated with JYPERK, instead of Cij, and do not 
take into account the backend. In combining CARMA 
data with single-dish data, we can overwrite JYPERK 
in CARMA data with dj , or define JYPERK for single- 
dish (NR045), as 



JYPERK 



2V2fcE 



'7?, CARMA 



(15) 



'7mb?7a^ V ^g,NR045 

Parameters are listed in Table [TJ 

5.3. Synthesized Beam Size 

The deconvolution process (e.g., CLEAN) replaces a 
synthesized beam with a convolution beam (typically a 
Gaussian). We determine the convolution beam size so 
that its beam solid angle matches that of the synthe- 
sized beam. Theoretically, the beam solid angle is an 
integration of a beam response function over 47r steradi- 
ans. In principle, we could calculate it by integrating a 
synthesized beam image or by taking the weight of the 
zero-spacing data (Appendix [B]). These methods worked 
reasonably well, but showed some error, introduced per- 
haps by the limited size of the beam image (not over 
47r steradians). In practice, we found that the following 
method provides better flux conservation: we calculate 
the total fluxes of the galaxy with the single-dish map 
and with the dirty image (with an unknown beam area 
as a free parameter) , and find the beam area that equal- 
izes these total fiuxes. The position angle and axis ratio 
of the beam is derived by a Gaussian fitting to the synthe- 
sized beam. The Gaussian is linearly scaled to reproduce 
the beam area from the flux comparison. 

If the solid angles do not match, the total fiux is not 
conserved in the final deconvolved map (e.g., CLEANed 
map). The CLEANed map has two emission compo- 
nents, deconvolved emission and noise/residual emission, 
and they have their own units of flux, J y/ (convolution 
beam) and Jy/ (synthesized beam), respectively. There- 
fore, the convolution beam smaller than the synthesized 
beam elevates the flux of the residual emission, while a 
larger beam reduces it. The error becomes particularly 
problematic for an object with extended, low-flux emis- 
sion (such as galaxies), which are inherently missed in the 
deconvolution process, but exist in the CLEANed map. 
The two units in the flnal map does not degrade an image 
quality much in case of the CARMA and NR045 image. 
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as long as the two beam areas are the same, because the 
synthesized beam is aheady similar to a Gaussian beam 
that we adopt as a convolution beam. 

We note that if the deconvolution procedure (such as 
CLEAN) can 'dig' all positive components down to the 
zero flux level, the convolution beam could have any 
shape. We also note that in case of pure interferome- 
ter observations, the beam solid angle is zero (Appendix 
[B|) . and thus, this method cannot be applied. 

5.4. Joint Imaging and Deconvolution 

The procedure for imaging and deconvolution is the 
same as the one in ^ 34. 31 but we add the term Wn/n in 
eq. ©, 

/(/, m) = W^h/h + Wolo + Wcic + Wn/n, (16) 

where /n and Wn are the image and weight from the 
NR045 visibilities, respectively. Wn is calculated with 
the pseudo primary beam Pn ( ^5.1|) and the theoretical 
noise ctn (= AS") derived with cq. (fTTI) . 

The number of NR045 visibilities is a control param- 
eter in this method (i i5.2p . and thus, should not be in- 
volved in weighting. Instead, the weight should be calcu- 
lated solely based on the sensitivity. We use the theoret- 
ical noise for weighting. If all visibilities have exactly the 
same fringe sensitivity, our weighting becomes the same 
as the conventional "natural" weighting. The sensitivity- 
based weighting for interferometer data was discussed in 
iBriggsl (|1995D . and our method extends it to the combi- 
nation with single-dish data. 

The robust weighting scheme su ppresses the pixels of 
high natural weights in uv space ffiriggsl ll995t ) - if the 
natural weight is lower than a threshold the weight is 
unchanged, but if it is higher, the weight of the pixel is 
set to this threshold. Our weighting scheme reproduces 
a natural weighting and works with the robust weighting 
scheme. We made two data cubes with robust = —2 
and -1-2. The resolution of the final combined data 
cube with robust = -2 is 3.7" x 2.9" (PA=79°) and 
5.08 kms-^ The RMS noise is 35 mJy/beam (i.e., 300 
mK) in lOkms^^ channel, robust = +2 gives the res- 
olution of 8.5" X 7.3" (PA=76°) and 5.08 km s'^, and 
the RMS noise of 52 mJy/beam (77 mK) in lOkms"^ 
channel. Figure [6] shows a synthesized beam pattern for 
robust — +2. Both cubes have the same total luminosity 
when the synthesized beam sizes are determined as in 

El 

6. INTEGRATED INTENSITY MAP AND IMAGE FIDELITY 

Figure [7] and [8] show the CO (J = 1-0) integrated 
intensity maps of M51, the combination of CARMA 
and NR045 data, with robust = -2 and -f2, respec- 
tively. The se maps are niade with the "masked moment 
method" in lAdler et al.l ()1992D . We also dropped the low 
sensitivity region (outer region) of the CARMA mosaic 
(see Figure [1]). The data with robust — -2 is used in 
following discussions, since it shows finer structures at a 
higher resolution. 

The combination of CARMA (15 antennas) and 
NR045 enables a full census of the population of giant 
molecular clouds (GMCs) over the entire galactic disk. 
Molecular gas emission in two spiral arms and interarm 
regions are prominent in this map. ,Koda et al.. (|2009i) 



showed the distribution of GMCs both in spiral arms and 
interarm regions, and the high molecular gas fraction in 
both regions. These two results suggest that stellar feed- 
back is inefficient to destroy GMCs and molecules, which 
is sup ported by a recent analysis by ISchinnerer et al.l 
(j2010t) . Molecular structures in the interarm regions were 
often an issue of de bate in previous observations due to 
poor image fidelity (iRand fc Kulkaf^ll990l: lAalto et al.l 
1999; Heifer et al.ll2003[ ). Figured compares the CO dis- 
tribution with a _B-band image from the Hubble Space 
Telescope (HST) and an 8/.*™ image from the Spitzer 
Space Telescope. Dust lanes in the i?-band image in- 
dicate the distribution of the dense interstellar medium 
(ISM), and the 8/zm image shows the distribution of the 
PAH (large molecules) illuminated by UV photons from 
surrounding young stars. The CO emission coincides 
very well with the dust lanes and 8/im emission in both 
spiral arms and interarm regions, which evidences the 
high image-fidelity over a wide range of flux. 

Figure [10] shows the NR045 map (left) and the ratio 
of the combined map (smoothed to ~ 20" resolution) 
over the NR045 map, i.e. recovered fiux map (right). 
The recovered fiux map shows an almost constant ratio 
^ 1 over the entire map, and no correlation with galac- 
tic structures (i.e., no size dependence - in contrast to 
the dependence expected in pure- interferometer maps). 
Some extended CO emission is not significantly detected 
at the high resolution of the combined image, but be- 
comes apparent when the image is smoothed. Note that 
the companion galaxy NGC 5195 is not included in the 
CARMA velocity coverage, nor in the combined map, 
though it is in the NR045 map. 

Both the main galaxy NGC 5194 and companion 
galaxy NGC 5195, are observed with NR045. The total 
flux of NGC 5194 is (1.022 ± 0.002) x lO'^ Jy • km/s in the 
N R045 map, which is consistent with the measurement 
of lHelfer et all (|200l . With the Galactic CO-to-Ha con- 
version factor Xco = 1.8xlO^'^cm^^[K-km/s]^^ and the 
distance of 8.2 Mpc, the total molecular gas mass in NGC 
5194is 4.9xlO^Mo. The Xpo similar to the Ga lactic one 
is found in M51 recently (jSchinnerer et al.]|2010[ ). The to- 
tal fiux of the combined cube is also 1.0 x 10''' Jy • km/s, 
consistent with the NR045-only measurement. The to- 
tal flux and mass of NGC 5195 is 162 ± 4 Jy • km/s and 
7.8 X 10 ''Mq, respectively. The errors are based on the 
RMS from the map, and do not include the systematic 
error due to the fiux calibration in the CARMA obser- 
vations (~ 15%). 

7. REQUIREMENTS 

There are requirements for sampling, field of view, uv- 
coverage, and sensitivity for single-dish data to be com- 
bined with interferometer data in an optima l manner. 

Fir st, a spatial fine sampling is necessary ()Vogel et al.l 
Il984f ). The half-beam sampling, a typical practice in 
most single-dish mapping observations, is not sufficient, 
since the aliasing effect destroys visibilities both at long 
and very short baselines. Figure [3] illustrates the ef- 
fect schematically: if the spatial sampling is 10.3" (= 
Acq / 52m, a typical sa mpling in NR045 observations; 
e.g. iKuno et al.l 120071) . the tail of the uv distribution 
leaks into baselines as short as ~ 7 m. Hence, the 
Nyquist sample of 11.9" (= Ago /45m) is necessary to 
properly reproduce visibilities up to the 45 m baseline. 
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The observing grid and pixel size in the NR045 map 
must be at most 5.96" (Figure |3]). 

The single-dish map should cover an area larger than 
the area of the joint map. The deconvolution with the 
single-dish beam (i 35.1.1|) causes artifacts at the edges of 
the images. It is ideal to have extra-margins with the 
width of a few single-dish beam sizes at each image edge. 

The sensitivity match between single-dish and inter- 
ferometer data should also be considered in matching 
their uv coverages; the maximum effective NR045 base- 
lines are limited by the matched sensitivity in our ob- 
servations. Only the baselines of about 1/4-1/3 of 
the 45m diameter take practical effect in the combina- 
tion. It is often discussed that a single-dish telescope 
needs to be about twice larger than the shortest baseline 
used in interferometer observations due to uncertainty 
in the single-dish beam shape and erro r s in pointing 
(see iKurono. Morita &: Kamazakil 120091 : iCorder et al.l 
120101 and references therein). In our case, the max- 
imum effective baseline is shorter than this length. In 
practice, interferometer data rarely cover the theoreti- 
cal minimum baseline (i.e. dish diameter). The long 
baselines of single-dish data do not have a sensitivity 
comparable to interferometer's (i i5.1|) . To avoid a gap 
in uv coverage without sensitivity loss, the diameter of 
the single-dish needs to be 3-4 times larger, unless the 
receiver of the single-dish telescope has a significantly 
higher sensitivity. The sensitivity match is discussed in 
detail in Appendix [C] 

8. COMPARISONS WITH OTHER METHODS 

Several methods for the combination of single-dish 
and interferometer data have been applied at millimeter 
wavelength. None of the previous data, however, have 
a sufficient overlap between single-dish and interferome- 
ter uv coverages (in the sense discussed in i 35.1.ip . The 
weighting schemes are artificial, rather than based on 
the sensitivity (i.e., data quality). Nevertheless, these 
methods have some advantages in simplicity, as well as 
d isadvantages in deta il. 

'St animirovic et ahl (|1999D introduce a combination 
method in the image domain. This method is adopted for 
the BIMA Survey Of Nearby Galaxies (BIMA-SONG) to 
comb ine BIMA int erferometer with NRA012 single-dish 
data ()Helfer et al.l l2003f). They set the weights to be in- 
versely proportional to the beam area (i.e., one term in 
eq. |9] and [TT\\ and add the dirty maps and beams of 
BIMA and NRA012 linearly (eq. [2]) to produce a joint 
dirty map and beam. The relative weights are manually 
and continuously changed with uv distance. The joint 
dirty map is then CLEANed with the joint synthesized 
beam. This method starts the combination process from 
images, rather than visibilities, and is simple. It should 
be able to use a more natural weighting scheme (e.g., 
sensitivity uw-distribution based on the beam shape and 
e g. see also App endix [C)) if software is developed. 

iWeifj et al. I ()200H ) also combine a single-dish map and 
CLEANed interferometer map. They deconvolve the 
single-dish map with its beam pattern and convolve the 
result with an interferometer convolution beam, so that 
the beam attenuation becomes the same for both single- 
dish and interferometer images. Then, they Fourier- 
transform both images and replace the interferometer 
data with the single-dish data at the central wti-spacing. 



CLEAN is performed separately for interferometer data 
alone, which does not take advantage of the high image 
fidelity of the combined map. Having only one control 
parameter - the choice of uv range to be replaced - can 
be advantageous. 

Visibilities are generated fr om a sing le -dish map 
by seve ral authors (Vogcl et al.l 119841: iTakakuwal 
200l iR odriguez-Fernandez. Petv fc Guetbl [20081 
Kurono. M orita & Kaniazaki 206W. Our method is in 
this branch. iHelfer et al.l (|2003) summarize difficul- 
ties to set the weights for this combination scheme, 
and conclude tha t it is too sensitive to the choice 
of pa ramete rs. IRodriguez-Fernandez. Petv fc GuethI 
(|2008D and IKurono. Morita fc Kamazakil (|2009() sug- 
gest to set the relative weight to obtain a cleaner 
synthesized beam shape, which is advantageous in 
deconvolution (e.g., CLEAN, MEM). More specifically, 
IRodriguez-Fernandez. P etv & Gucth, (|2008[ ) set the 
single-dish weight density in uv space equal to that 
of the interferometer vi sibilities that surround the 
single -dish uv coverage. IKurono. Morita fc Kamazakil 
(|2009| ) adjusted the relative weight to zero out the total 
amplitude of the sidelobes of a synthesized beam. Our 
weighting scheme is more intrinsic to each set of data 
and is based solely on their qualities; the single-dish 
weight is independent of the interferometer data and is 
set based on the RMS noise of a single-dish map. The 
weight is not a parameter of choice. 

In pure interferometer imaging, the synthesized beam 
shape is historically controlled by changin g the weight 
density in uv space. The robust parameter ( Briggsiri995fl 
is a famous example that converts the weight smoothly 
from the natural to uniform weightings. Once the weight 
is set in our method, the robust weighting works even for 
the combined data, exactly as designed for pure interfer- 
ometer data. 

9. SUMMARY 

We describe the CARMA observations at the early 
phase of its operation, and the OTF observations with 
the multi-beam receiver BEARS at NR045. The stan- 
dard reduction of CARMA and NR045 data are also 
discussed and extended to the combined data set case. 

We explain the basics of the imaging technique for het- 
erogeneous array data, and show that the combination 
of interferometer and single-dish data is an extension of 
the imaging of heterogeneous array data. We introduce a 
method of combination of interferometer and single-dish 
data in uu-space. The single-dish map is converted to 
visibilities in uu-space. The weights of the single-dish 
visibilities are determined based on the RMS noise of 
the map, which is more natural than any other artificial 
weighting schemes. The synthesized beam size is deter- 
mined to conserve the flux between the dirty beam and 
the convolution beam. Comparisons with other meth- 
ods are discussed. The advantages and disadvantages of 
those methods are summarized. In the appendices, we 
discuss the matching of single-dish and interferometer 
sensitivities for the combination of the data. 

The resultant map shows the high image fidelity and 
reveals, for the first time, small structures, such as giant 
molecular clouds, both in bright sp iral arms and in faint 
inter- arm regions (Koda et al.*'2009^. From the new map, 
we calculate that the total masses of NGC 5194 and 5195 
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are 4.9 x 1O^M0 and 7.8 x 1O^M0, respectively, assuming 
Xco = 1-8 X 102"cm-2[K-kni/s]-i. 

The combination method is designed on a platform of 
available software (i.e., MIRIAD) and generates a finite 
number of discrete visibilities from a single-dish map. Fu- 
ture software should enable data manipulation directly 
on maps (grids) both in real and Fourier spaces, instead 
of in visibilities f Appendix iPj). The weights can be deter- 
mined on a grid basis, rather than on a visibility basis. 
Even in such cases, the weights should be determined 
from the RMS noise of the map which are related to the 
quality of data. 
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Figure 1. Observed areas of M51 over the Spitzer 24/jm image of M51. Red crosses: 151 pointing positions of CARMA observations. 
Blue contours: sensitivities (RMS noise) of 120, 140, 160, 180, and 200 % v^ith respect to the central maximum sensitivity. Most part over 
the galactic disk has a uniform sensitivity. Black box: NR045 coverage. 
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Figure 2. The 

CARMA (black) 



uv coverage at the central pointing in unit of kilo-lambda. Left: CARMA uv coverage. Right: the central region of 
and NR045 (red) uv coverages. 
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Figure 3. Schematic illustrations of the aliasing effect in the uv coverage of the NR045 map. The alias of the NR045 data in uv space 
destroys not only long, but also short baselines. In the case of 10.3" sampling, the baselines as short as 7 m are affected by the alias. 
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Figure 4. Flow chart of the process of combination. The NR045 cube is deconvolved with the NR045 beam, multiplied with an arbitrarily 
selected primary beam, and Fourier-transformed into Mti-space. The transformed data are re-sampled with a Gaussian uv distribution to 
produce NR045 visibilities. The weight of the NR045 visibilities, with respect to CARMA, are determined based on the RMS noise of the 
NR045 cube. 
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Figure 5. Sensitivities as a function of baseline length for NR045 (diamonds) and CARMA (crosses). The bottom panel is the same as 
the top panel, but only for short baselines. Vertical lines, from right to left, correspond to the antenna diameter of NR045 (45m) and its 
three quarter, half, and quarter lengths. 
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Figure 6. Synthesized beam, i.e., the combination of the 4 baseline types (CARMA, HATCREEK, OVRO, and NR045), for robust=+2. 
Left: Synthesized beam pattern. Contours are -4, -2, 2, 4, 10, 20, 40, 60, 80, 100 % of the peaJf. Dashed lines are negative contours. 
Middle and right: Slices of synthesized beams along RA and DEC directions, i.e. a superposition of all the beams at the 151 pointings. 
The synthesized beam patterns, i.e., uv coverages, are very similar for all the pointings. 
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C0(,/ = 1-0) Imaging of M51 with CARMA and NR045 
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Figure 9. CO contours on an HST i?-band image (left) and on a Spitzer 8fim image (right). The CO image is first smootlied to a 30" 
resolution to mask out < 3cr pixels in the smoothed image. The CO contours are made with the masked robust=-2 map (not with the 
integrated intensity map), and are the superposition of the 1.5, 5.0, 10.0, 15.0 a contours in all lOkms"^ channels. The extended emissions 
at the 1.5cr level are significant in the smoothed image. The HST image is derived by dividing the original image by an axisymmetric 
luminosity profile to visualize the dark dust lanes. The CO emission coincides with the dust lanes both on bright spiral arms and in dimmer 
interarm regions, indicating the high fidelity of the CO data. The Spitzer Sfim image traces the dense interstellar medium illuminated by 
UV photons, and the CO coincides with the Sfim emission as well. The circle at the lower-right corner of the left panel has the diameter 
of 4" (roughly the size of beam) . 
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Figure 10. Left: NR045 CO ( J = 1 — 0) map. Right: Recovered flux map, i.e., the ratio of the combined map over NR045 map. The 
combined map {robust = —2) is smoothed to ~ 20" resolution for the comparison. The ratio is 1 over the entire map, and the flux 
recovery is very good. The companion galaxy NGC 5195 is not in the CARMA velocity coverage. Note that the error in the ratio map 
varies across the map, depending on the brightness of emission. 
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Figure 11. Sensitivity distribution with baseline length for CARMA (solid lines) and NR045 (dashed lines). The labels are the imaging 
sensitivities of CARMA and NR045. The minimum and maximum baseline lengths of CARMA are 10 and 300 m, respectively (for C & 
D configurations). The sensitivities, pixel sensitivity AS^ (y-axis) and imaging sensitivity A5iinage(= ^'S''), are calculated for the velocity 
channel width of A^; = lOkms"'^. 
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Table 1 

Antenna Parameters at 115 GHz 



Parameter Unit Hatcreek (6m) OVRO (10m) CARMA (6m-10m) NR045 



Main Beam Size FWHM arcsec 100 60 77.5" 19.7* 

Beam Solid Angle arcsec^ 1.51 X lO" 6.77 X 10=* 1.01 X 10"'" 1.10 X 10=*' 

Quantum Efficiency 77q ■ ■ ■ 0.87 0.87 0.87 0.87 

Main Beam Efficiency ry^b ■ ■ ■ 0.41 0.61 0.50" 0.40 

Noise Coef. (general) Cij Jy/K 116.8 52.2 78.1" 12.0 

Noise Coef. (MIRIAD) JYPERK Jy/K 145.3 65.0 97.2" 14.9 



Note. — Table [T] 

Geometric mean of Hat creek and OVRO values. 
After re-gridding ( ^3^ . 
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APPENDIX 
WEIGHT FUNCTIONS 

The dirty image /"^"^ and synthesized beam B are defined with a set of visibihties V{u,v) as 

I'^"\l,m) = J J V{u,v)W{u,v)e^^''-''^+^"'Uudv (Al) 

and 

B{l,m)= [ /w^(M,v)e2"("'+^'")dudz;, (A2) 



where (u, v) is the coordinates in the uv space. 
The samphng and weighting function of visibihties W{u, v) can be written more exphcitly as 

M 

W{u,v) ^Y.RkTkDkS{u - Uk,v - Vk), (A3) 

where Tk is the tapering function, and Dk is the density weighting function (see iTavlor. Carihi fc Perlevl[T999l) . M is 
the number of visibihties obtained in observations. Tk and Dk are arbitrary functions, and are often used to control 
the synthesized beam shape and noise level. For example, the Gaussian taper is Tk = exp(— -y/u^, + v'f, '/2a^) with the 
half power beam width 6'hpbw = ■\/2 ln2/7r/a ~ 0.37 /a [radian]. The natural and uniform weightings are Dk — 1 and 
Dk = i/Nk, respectively, where Nk is the number of visibilities within a pixel in uv space. 

Rk is a weight based on noise, and has the relation Rk = l/AS^ with ASk (= AS-^ in §4.4.ip . The theoretical noise 
of an image a can be calculated as 



AI \ / M \ M 



E ^'DiRk E / E (A4) 



BEAM SOLID ANGLE 
The beam solid angle D,a of a synthesized beam (dirty beam) is defined as 



VtA= B{l,m)dldm (Bl) 



W{u,v) 



dudv (B2) 



^W{0,0) (B3) 

Eq. (IA2[) is used bewteen eq. (|B1I) and (jB2[) . The bracket in eq. (|B2p is a (5- function. We assumed that the maximum 
of B{1, m) is normalized to 1. 

Pure interferometer observations do not have zero-spacing data, and therefore, JIa = 0. A Gaussian beam B{1, m) = 
exp[— (P -|-m^)/2CT^] has W{u,v) = 27rfT^ exp[— (m^ + v^)/2ap], where crp = l/^na. Therefore, VIa = 2na'^. 

SENSITIVITY MATCHING BETWEEN CARMA AND NR045 

Matching the sensitivities of CARMA and NR045 is crucial in combination. The sensitivity requirements of the 
interferometer and single-dish maps are important for the observation plan, and a simple way to calculate matching 
sensitivities is therefore important. 

One approach is to match the sensitivities in uu-space around the uv range (baseline range) where two data sets 
overlap. In other words, we want to match the pixel sensitivities AS^ of CARMA and NR045, i.e, their sensitivities 
at pixel (u, v). Among some definitions of sensitivity (e.g., !j4.4.1l) . the imaging sensitivity AS^, i.e., noise fluctuation 
in the map (eq. [9]), is most used to characterize the quality of map and to estimate the feasibility of observations. 
Therefore, we first derive the relation between the imaging and pixel sensitivities in uv space. 

For simplicity, we assume that all CARMA antennas are identical to each other, having exactly the same Tgys and 
Cij . Then, a sensitivity is 

A5.C„-^, ,C1) 

where B is the channel width and t is the integration time. is applied to both CARMA and NR045, and can 
mean one of the following three sensitivities: fringe sensitivity AS*^ when t is the integration time of a visibility tvis 
(eq. [7]); imaging sensitivity A5" when t is the total integration time, i.e, tvisNyis, where TVvis is the total number of 
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visibilities; and pixel sensitivity AS'p when t is the total integration time of the pixel at {u,v), i.e, tvisn{u,v), where 
n{u, v) is the number of visibilities in the pixel. Therefore, the imaging and pixel sensitivities are related as 



AS^u,v)=AS\ (C2) 
y n(u, V) 

Hereafter, we derive the relation between A'vis and n{u, v). 

The n{u, v) for the interferometer (e.g., CARMA) was discussed bv lKurono. Morita fc Kamazakil (120091 ) . For a target 
at a reasonably high declination, synthesis interferometric observations provide a visibility distribution of n{b) cx 1/6, 
where b is the uv distance b — y/u'^ + v'^. Th visibilities (total of A'vis) are distributed within the minimum and 
maximum baseline lengths, by^in and 6max respectively. From A'vis — n{b)2iibdh, we derive 

n{b) = — (C3) 

27r(0,„ax — Omin) 

The n{u, v) for a single-dish telescope (e.g., NR045) is determined by the beam shape of the telescope. We assume 
a Gaussian beam shape, cx exp[— (P -f TO^)/2cr^], in sky coordinate {l,m). The full width half maximum (FWHM) 
of the beam is FWHM = 2\/21n2cr. The n{u, v) is proportional to the Fourier transformation of the beam shape, 
cx exp [-(27rcr)262/2]. The total of A'vis visibilities are within the uv range from zero to the antenna diameter d. Thus, 

"^^)= l-exp[-(2..W2] ^ ' ■ (^'^ 

Eq. (jC2p . (jC3p . and (jC4p give the pixel sensitivity for interferometer A5'?Jj^(6) and single-dish A5|'^(6). Equalizing 
these two ASl^^{b) — AS^^{b) at b = ^overlap where the two uv coverages overlap leads to a relation between the 
image sensitivities (i.e., RMS map noise) of the interferometer and single-dish. This relation is a rough measure of 
the matched sensitivities for the combination of the single-dish and interferometer, and would be useful in planning 
observations. The sensitivity matching can be calculated more accurately with eq. (jC2p . as performed in i jS.l.ll if we 
know an accurate uv coverage n{u, v) of interferometer observations. 

Figure [TT] plots the pixel sensitivities AS^ of CARMA and NR045 as function of baseline length b for fixed image 
sensitivities AS*'. We set 6,nin and b^^x to 10 and 300 m (~ 4 and 115kA), respectively, for CARMA C & D- 
configuratlions. The CARMA and NR045 uv coverages overlap significantly between 4 and 10 kA. The NR045 
noise (sensitivity) increases rapidly beyond the baseline length of about half the diameter (~ 8kA) , and CARMA can 
complement the uv range beyond that. The imaging sensitivities of our CARMA and NR045 observations are 27 and 
155 mjy in the velocity width of 10kms~^, respectively. Their sensitivities match around b ^ 4-6kA, within the range 
where the uv coverages overlap. 

APPLICATION TO GRID-BASED COMBINATION SCHEME 

The new combination technique discussed in this paper converts a single-dish map to a finite number of visibilities 
(discrete data points in uw-space). The weight of each single-dish visibility is determined based on the RMS noise of 
the map (i.e., the quality of the data) using the fringe sensitivity AS^ (eq. I13[) . In the Fourier transformation, the 
visibilities are mapped to a grid in uu-space, and the pixel sensitivity AS'P is calculated for each pixel of the grid by 
summing up the AS*^ of all the visibilities in the pixel. The AS^ for single-dish data can be calculated directly without 
going through visibilities once proper software is developed. 

The pixel sensitivity AS'P for the single-dish data can be defined with eqs. (jC2p(|C4p . The RMS noise of the single- 
dish map AS' gives the normalization of the equations. Eq. (jC4p is for a Gaussian beam, and could be replaced with 
some other shapes, such as a Fourier transformation of a single-dish beam or PSF if we have better knowledge of them. 
The ASP for the interferometer data should be calculated from the fringe sensitivities of visibilities AS^ using Cij (eqs 
[3 [8]) - mapping the visibilities onto a grid in uv-space and summing up the fringe sensitivities in each pixel. These 
ASP naturally set the relative weight of the single-dish and interferometer data. 



